Distribution of the local density of states, reflection coefficient and Wigner delay time 
in absorbing ergodic systems at the point of chiral symmetry 
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Employing the chiral Unitary Ensemble of random matrices we calculate the probability distribu- 
tion of the local density of states for zero-dimensional ("quantum chaotic") two-sublattice systems 
at the point of chiral symmetry E — and in the presence of uniform absorption. The obtained 
result can be used to find the distributions of the reflection coefficent and of the Wigner time delay 
for such systems. 
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Spectral and transport properties of quantum systems 
with chiral symmetries attracted a considerable attention 
recently Q, [111 EH HUH. The systems of the discussed 
type are characterized by random Hamiltonians of the 
form: 



H 
fft 



(1) 



serving to describe various physical situations, e.g. 
lattice systems with bond disorder of pure nearest- 
neighbour hopping typeQ, 0,0,0], random magnetic flux 
models appearing in the context of the Quantum Hall 
effect^], as well as lattice QCD models @|. Many other 
interesting examples can be found in the cited literature. 

Replacing the block H in Eq.Q with random rectan- 
gular matrices of the size Na x Nb with independent, 
Gaussian-distributed complex entries we arrive at the 
object known as the chiral Gaussian Unitary Ensemble 
(chGUE) 0] . In the language of the two-sublattice disor- 
dered lattice system this model is known to correspond 
to the so-called universal ergodic (or " zero-dimensional" , 
or " quantum dot" ) limit adequate if we are interested in 
energies smaller than the so-called Thouless energy E c . 
The latter requirement means that the time to diffusively 
propagate through the system is much shorter than %/ A, 
with A standing for the mean level spacing in the rele- 
vant energy range. Spectral statistics of chGUE matrices 
is well-studied by various methods and many explicit re- 
sults are available EEEEEGjI- 

One of the most striking features of the spectra of chi- 
ral systems is a specific behaviour of the spectral char- 
acteristics at zero energy. In zero-dimensional chiral sys- 
tems this specific behaviour is manifested, in particular, 
via vanishing of the mean eigenvalue density p av (E) at 
the middle of the band, E = 0, as long as Na — Nb- 
For Na ^ Nb there exist exactly \Na — Nb\ zero eigen- 
values. All nonzero eigenvalues for chiral models appear 
in pairs ±A„, with the corresponding eigenvectors being 

! Here u„ and v„ are eigenvectors of the 



and Nb X Nb, respectively. To ensure the orthonormal- 
ity of the eigenvectors (w^w^) — S mn we choose the 
normalisation vJjV„ = u£u n = I/a/2- 

We first assume, for simplicity, that the number of lat- 
tice sites in the two sublattices are equal: Na = Nb = N 
and define in a usual way the exact local density of states 
(LDOS) at a site ja — 1, 2, N of the sublattice A via 
the relation: 



Pj A ( E v) = -Mi-al „ 1 n, 
IT E q — n 

9 N W 
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where u n (j) stands for j — th component of the eigenvec- 
tor u n . Here we assume that the energy E n = E — if] 
is complex, the parameter rj > standing for the level 
broadening due to an absorption, i.e. uniform-in-space 
losses of the flux of particles in the sample. In particu- 
lar, for the energy E chosen precisely in the center of the 
spectrum ("the point of chiral symmetry") the LDOS is 
given by the most simple expression: 



>1 



(4) 



For the lattice sites of the sublattice B the LDOS is given 
by the same formula, with u„(j) replaced by v n (j). 

The above expression is a very convenient starting 
point for calculating the distribution function of the 
local DOS at the spectrum center. In doing this we 
employ the method suggested for GUE without chiral 
structure in|lfij]. Namely, we consider the generating 
function (the Laplace transform of the probability den- 



sity) P{t) 



-hGUE wnere the angular brack- 



Hermitian matrices HW and WH, of the size Na X Na 



ets stand for the ensemble averaging. Using that the 
eigenvalues A„ are statistically independent of the eigen- 
vectors u„,v„, and the components of the latter vec- 
tors behave in the large N limit like independent com- 
plex Gaussian variables with zero mean and variances 
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{\ u n{j)\ 2 ) ch cu e = one can ^ rs * P er f° rm the ensem- 
ble average over eigenvectors. After simple manipula- 
tions one arrives at the following representation for the 
generating function: 



P(t) 



' det (V + H H 
, dct (fj 2 + HH^ 



(5) 



H 



where we denoted fj = \Jr\ 2 + i]t/(irN). In this way the 
problem of calculating the generating function in ques- 
tion is reduced to evaluating the ensemble average of the 
ratio of characteristic polynomials of the random matri- 
ces HW . The latter object (and its generalisations) were 
intensively studied in the literature recently, see |13t llfij , 
and the result is readily available: 

P(t) = \ [eh(e)K (Ve(e + 2i/7r)) (6) 

+ (Ve(e + 2*/tt)) Ki (V^+W^)) 'o(e) 

where we denoted e — 2Nr/ and considered this parame- 
ter to be finite when performing the large-N limit. Here 
I n (z),K n (z) stand for the modified Bessel and Macdon- 
ald functions of the order n, respectively. The distribu- 
tion function V(pa) of the LDOS can be found by invert- 
ing the Laplace transform Eq. JBJ) with respect to t and is 
given by: 



2pA 



h(e) + — J (e) 
irpA 



(7) 



Having the LDOS distribution at our disposal, we can 
further utilise it for obtaining distributions of varous 
physical quantities of much interest describing scattering 
of a quantum particle from the two-sublattice disordered 
sample. In doing this we follow the method suggested 
recently by the first author in [l^. Consider a particle 
injected into the site ja of the sublattice A through a 
one-channel infinite lead perfectly coupled to the system. 
The scattering of such a particle from our absorptive sam- 
ple is described by the scattering matrix SjA expressed in 
terms of the matrix elements of the resolvent associated 
with the chiral Hamiltonian TL: 



S jA {E n ) - - 



1-iG 



iG 



G = 7T(j A \ 



1 



En —H 



A ■ 



(8) 



In this way the problem of investigating statistics of the 
scattering matrix amounts to that of the diagonal entry 
G of the resolvent. For an arbitrary energy E the vari- 
able G is a complex quantity and one has therefore to 
know the joint probability density of its real and imagi- 
nary parts. However, at the point E — of chiral sym- 
metry G turns out to be purely imaginary and therefore 
just proportional to the LDOS, Eq.(@J, whose distribu- 
tion T(pa) was calculated above, Eq.Q). As a result we 
can calculate the distribution of the reflection coefficient 



p ,(Pa) 
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FIG. 1: Numerical (symbols) and analytical (solid lines) re- 
sults for the LDOS distribution of the sublattice A for e = 1, 
N B = 300 and N A = 300-1, 1 = 0, 1, 2, 3. 



P^jAiv) = \SjA(v)\ 2 - Below we present the distribution 
of a related object r = 1 — RjA — 47t/9ja/(1 + ^PjA) 2 
which, following [J^i we suggest to call the "probability 
of no return" (PNR). This quantity has meaning of the 
quantum-mechanical probability for a particle entering 
the system via a given channel never exit back to the 
same channel. It is well-defined for a given realization of 
disorder and will therefore show sample-to-sample fluctu- 
ations. Let us stress that in the case rj = of no internal 
dissipation S— matrix is unitary for any energy and PNR 
is vanishing identically. For a system with absorption the 
unitarity is violated, giving rise to a nontrivial statistics 
of t: 



V(t) 



-VT 



7i(e)+(--l 



(9) 



The above distribution can be further utilized for 
finding statistics of such interesting characteristic of 
quantum scattering as the so-called Wigner delay time 
Tw = —idS/dE intensively studied in recent years, see 
, , [H , H3 and references therein. This can be done 
again following |17|. see also |2lJ. Indeed, remembering 
that nonzero values of r arise solely due to an absorp- 
tion, one can relate for small absorption the PNR value 
to the Wigner delay time by r ~ 2ervK0- Taking the 
corresponding limit in the distribution J5J| we arrive at 
the distribution of the scaled variable tw = tw/2N: 

1 l_ 

V(?w) = — e ^ (10) 
T w 

The expression Ea. (|10|l ) should be compared with the 
corresponding formula for the one perfect channel at- 
tached to an ergodic system modelled by GUE, i.e. with- 
out the chiral symmetry 0: V(fw) — -4-e The 
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FIG. 2: Numerical (symbols) and analytical (solid lines) re- 
sults for the LDOS distribution of the sublattice B for 6 = 1, 
N B = 300 and N A = 300 - I, 1 = 0, 1, 2, 3. 



FIG. 3: Numerical (o) and analytical (solid lines) results for 
the distribution of scaled Wigner delay times for sublattice 
A, for e = 1, N B = 300 and N A = 300 — Z, 1 = 0, 1, 2, 3. 



crossover between the two types of behavior is expected 
to happen in the present model with increasing energy 
parameter E from zero to values much exceeding the 
mean level spacing. Investigating this interesting ques- 
tion requires however much more elaborate methods and 
is left for future research. 

Let us now briefly discuss generalization of our results 
to the case of non-equivalent sublattices N B > Na- Note 
that in this case exactly Nb — Na — I zero modes are 
supported by the sublattice B, whereas only modes with 
nonzero eigenvalues give contribution to the LDOS for 
the sublattice A. Assuming that I <C Na,Nb one finds 
that LDOS distribution for the sites belonging to the 
sublattice A is given by formula generalizing Eq.Q: 



Vi(Pa) 



iJ+i 



-%(npA+l/irpA) 



1 



KPA 



Ii(e) 



(11) 

whereas the LDOS distribution Vi(pb) for sites belonging 
to the sublattice B is obtained from the above formula 
by replacing I with —I. In Fig. 0) and Fig. (J2J) we com- 
pare this prediction with the results of direct numerical 
simulations of the random matrix Hamiltonian Eq.Q), 
for sublattices A and B correspondingly. It is instructive 
to compare the behavior of the mean LDOS for the two 
sublattices in the limit e — » 0: 

(Pa) = | [ii+itoJS-xCe) +Ii(e)K l (e)} -> ± (12) 

(pb) = \ [I,-i(e)tf I+1 (e) +Ii(e)K l (e)] - ~ (13) 

The singularity in the latter expression reflects presence 
of exactly I zero modes supported by the sites of sublat- 
tice B. 



After straightforward manipulations one can find that 
the distribution Eq.(|5J) is replaced by the following two 
formulas: 



Vl(r A ) = 



T A \J\ - T A 

x [(< + ^) We) + K +1 + ^)] 



Vi{t b ) = 



[(v l + + v[ 



n-i 



T B \J\ - TB 

(e) + (v l + - 1 +v l S 1 )l l (e)]. 



(14) 



(15) 



valid for the sublattices A and B, respectively. Here the 
parameters v+,V- satisfy the conditions: v + v- = 1 and 

i(u++u_) = ( T J^ B) — l) • Finally, performing in Ea. H14|) 
the limiting procedure discussed above, we arrive at the 
distribution of the (scaled) Wigner delay time fw for the 
one-channel lead attached to a site belonging to the A— 
sublattice: 



Vi(t w ) 



1 1 



(16) 



WA 



which generalizes the expression Ea. (|10|) . In Fig. (J3j) we 
present the distribution of the Wigner delay times calcu- 
lated numerically for various values of I compared with 
this analytical result. 

In fact, the time-delay distributions Eqs. (|10fl and i|16|) 
can be obtained in a more direct way. For this we no- 
tice that the expression Eq. JSJ) for the unitary scattering 
matrix Sa(E) at real energies E can be written in an 
equivalent form: 



S A {E) 



det E 2 - iEP. 



3A 



det ( E 2 + iEP, 



(17) 
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where we used the Na x Na projection matrix Pj A — 
\]a >< ]a\- Calculating now the Wigner delay time as 
the energy derivative of the logarithm of the above equa- 
tion, and setting E = we arrive after straightforward 
manipulations to a simple formula: 

TW =2<j A \(H&y 1 \3 A >=2Y:^^ (18) 

n=l n 

Comparison of this expression with Eq. (0J makes it clear 
that the Laplace transform P(s) of the probability den- 
sity for t\y — tw/2Na can be obtained via the same 
manipulations as those leading to Eqs.© and ©• For 
the general case < / = Nb — Na <C Na, Nb we find: 

/ det (hhA \ 2 
\det [s/Nl + HH^f U 

which can be immediately inverted yielding the distribu- 
tion Eq.USJl. 

It is easy to repeat the same manipulations for the lead 
attached to a site of the larger sublattice B and to find 
that Eq. (|18|l is replaced formally with 

tw = 2<5b\(&h) X \3b > (19) 

This expression, however, does not make any sense be- 
yond the case I = of equivalent sublattices, since the 
corresponding Nb x Nb matrix WH is singular. We 
conclude therefore that the Wigner delay time at zero 
energy diverges for the sites belonging to the larger sub- 
lattice due to presence of zero modes. 

In summary, we presented analytical as well as nu- 
merical results for the distribution of the local density of 
states for zero-dimensional chaotic system at the point of 
chiral symmetry E = 0, and further utilized it for study- 
ing statistics of the reflection coefficient and Wigner delay 
times. Apart from investigating a gradual breakdown of 
the chiral symmetry with increasing energy parameter E 
in ergodic systems, an interesting and challenging ques- 
tion is to study the LDOS fluctuations in chiral systems of 
higher dimensions. In particular, in disordered quantum 
wires the mean DOS is known to display a logarithmic 
singularity 0, 0- It would be interesting to understand 
how that singularity affects the LDOS distribution, and 
how it is reflected in fluctuations of the related scattering 
characteristics. 
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